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Motivated by the recent experimental discoveries of multi-quark candidates, e.g., 
■ the ©+(1540), we study multi-quark systems in lattice QCD. First, we perform 

' accurate mass measurements of low-lying 5Q states with J = 1/2 and I = in 

' both positive- and negative-parity channels in anisotropic lattice QCD. The low- 

est positive-parity 5Q state is found to have a large mass of about 2.24GeV after 
the chiral extrapolation. To single out the compact 5Q state from NK scatter- 
ing states, we develop a new method with the hybrid-boundary condition (HBC), 
and find no evidence of the compact 5Q state below 1.75GeV in the negative- 
parity channel. Second, we perform the first study of the multi-quark poten- 
tial in lattice QCD to clarify the inter-quark interaction in multi-quark systems. 
The 5Q potential V5Q for the QQ-Q-QQ system is found to be well described by 
the "OGE Coulomb plus multi-Y Ansatz": the sum of the one-gluon-exchange 
(OGE) Coulomb term and the multi-Y-type linear term based on the flux-tube 
picture. The 4Q potential V4Q for the QQ-QQ system is also described by the 
OGE Coulomb plus multi-Y Ansatz, when QQ and QQ are well separated. The 
4Q system is described as a "two-meson" state with disconnected flux tubes, when 
the nearest quark and antiquark pair is spatially close. We observe a lattice-QCD 
evidence for the "flip-flop", i.e., the flux-tube recombination between the connected 
4Q state and the "two-meson" state. On the confinement mechanism, the lattice 
QCD results indicate the flux-tube- type linear confinement in multi-quark hadrons. 



"Talk presented by Fumiko Okiharu at "Quark Nuclear Physics 2005" (QNP05). 
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1. Introduction 

Recently, several new particles were experimentally discovered as the candi- 
dates of multi-quark hadrons. The + (154O) was found in 2002 at SPring- 
8 1 as a pentaquark (5Q) baryon and was confirmed in other experimental 
groups (ITEP, JLab, ELSA). 2 - 3 ' 4 The 9+(1540) has the baryon number 
B = 1 and the strangeness S = +1, and hence it is a manifestly exotic 
baryon and is considered to be a pentaquark (uudds) in the valence-quark 
picture. For other candidates of pentaquarks, the S~~(1862) (ddssw) was 
found at CERN, 5 and the C (3O99) (uuddc) was reported as a charmed 
pentaquark at HERA. 6 Subsequently, the candidates of tetraquark (4Q) 
mesons were observed. The A(3872) 7 ' 8 ' 9,10 was found in the process of 
B+ -^K+ + A(3872) K+ + tt-tx+J/^ at KEK. 7 The A(3872) is much 
heavier than the J/tp, and its mass is close to the threshold of D° (cu) 
and D°* (uc). Nevertheless, its decay width is very narrow as T < 2.3MeV 
(90% C.L.). These features indicate the A(3872) to be a tetraquark, e.g., 
a bound state of D° and D°*. Similarly, the D s (2317) n ' 12 is expected to 
be a tetraquark candidate. A common remarkable feature of all these ex- 
otic hadrons is their extremely narrow width, which seems a new mystery 
appearing in hadron physics. In any case, these discoveries of multi-quark 
hadrons are expected to reveal new aspects of hadron physics. 

In the theoretical side, the quark model is one of the most popular 
models to describe hadrons. In the quark model, mesons and baryons are 
usually described as qq and 3q composite particles, respectively. In more 
microscopic viewpoint, quantum chromodynamics (QCD) is the fundamen- 
tal theory to describe the strong interaction. In terms of QCD, not only 
ordinary qq mesons and 3q baryons, but also exotic hadrons, such as multi- 
quark hadrons (qqqq, qqqqq, • • • ), hybrid mesons (qqg, • • • ), hybrid baryons 
(<Z<Z<Z.9j ' ' ' ) an( i glueballs {gg,ggg, •••) are expected to appear. We here 
aim to study these multi-quark hadrons directly based on QCD. Even at 
present, however, it is rather difficult to deal with the low-energy region 
analytically in QCD owing to its strong-coupling nature. As an alternative 
way, the lattice QCD Monte Carlo simulation is established as the power- 
ful method to treat non-perturbative nature of low-lying hadrons including 
exotic hadrons. In this paper, we perform the following two lattice QCD 
studies to clarify the properties of multi-quark systems. 

First, we investigate the mass and the parity of the 5Q system in lattice 
QCD. As for the parity assignment of the 0(1540), it is not determined 
experimentally yet, and little agreement is achieved even in the theoretical 
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side. On one hand, the positive-parity assignment is supported by the chiral 
soliton model 13 and the diquark model. 14 On the other hand, the negative- 
parity assignment is supported by the nonrelativistic quark model, 15 the 
QCD sum rule 1 and so on. These predictions have been done with model 
calculations, but these models were originally constructed only for ordinary 
hadrons. In fact, it is nontrivial that these models can describe the multi- 
quark system beyond the ordinary hadrons. To get solid information for 
the multi-quark systems, we study their properties directly from QCD by 
the lattice QCD simulation, 17,18 which is the first-principle calculation and 
model independent. 

Second, we study the inter-quark interaction in multi-quark systems in 
lattice QCD. The inter-quark force is one of the most important elementary 
quantities in hadron physics. Nevertheless, for instance, no body knows the 
exact form of the confinement force in the multi-quark systems directly from 
QCD. In fact, some hypothetical forms of the inter-quark potential have 
been used in almost all quark model calculations so far. Then, the lattice 
QCD study of the inter-quark interaction is quite desired for the study of 
the multi-quark systems. It presents the proper Hamiltonian in multi-quark 
systems and leads to a guideline to construct the QCD-based quark model. 
In this paper, to clarify the inter-quark force in the multi-quark system, 
we study the static multi-quark potential systematically in lattice QCD 
using the multi-quark Wilson loop. We investigate the three-quark (3Q) 
potential, 19,20,21 which is responsible to baryon properties, and perform the 
first lattice-QCD study for the multi-quark potential, the tetraquark (4Q) 
and the pentaquark (5Q) potentials. 22 ' 23 ' 24 ' 25 ' 26 

We show in Fig. 1 our global strategy to understand the hadron proper- 
ties from QCD. One way is the direct lattice QCD calculations for the low- 
lying hadron masses and simple hadron matrix elements, although the wave 
function is unknown and the practically calculable quantities are severely 
limited. The other way is to construct the quark model from QCD. From 
the analysis of the inter-quark forces in lattice QCD, we extract the quark- 
model Hamiltonian. Through the quark model calculation, one can obtain 
the quark wave-function of hadrons and more complicated properties of 
hadrons including properties of excited hadrons. 

This paper is organized as follows. In Section 2, we perform accurate 
mass calculations for low-lying 5Q systems in anisotropic lattice QCD. 27 In 
Section 3, we perform the systematic study of the inter-quark interaction 
in multi-quark systems. 22 ' 23 ' 24 ' 25 ' 26 Section 4 is devoted for the summary 
and concluding remarks. 
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Figure 1. Our global strategy to understand the hadron properties from QCD. One 
way is the direct lattice QCD calculations for the low-lying hadron masses and simple 
hadron matrix elements, although the wave function is unknown and the practically 
calculable quantities are severely limited. The other way is to construct the quark model 
from QCD. From the analysis of the inter-quark forces in lattice QCD, we extract the 
quark-model Hamiltonian. Through the quark model calculation, one can obtain the 
quark wave-function of hadrons and more complicated properties of hadrons including 
properties of excited hadrons. 



2. Lattice QCD Study for the ©+(1540) 

There are many theoretical studies of the + (154O) in order to clarify 
its physical properties. 28,29 For instance, as was briefly mentioned in the 
previous section, its parity is not yet determined experimentally, and is 
still unsettled also in the theoretical side. Furthermore, the existence of 
the O + (1540) itself is not experimentally established as a pentaquark res- 
onance. In fact, recent some high-energy experimental groups reported 
no evidence of the + (154O), 30 ' 31 ' 32 while several low-energy experimental 
groups reported its existence. 

Also in lattice QCD, there is no consensus on the existence and the par- 
ity assignment of the + (154O). Two early works supported the negative- 
parity state for the + (1540), 33,34 while one early work supported the 
positive-parity state. 35 (To be strict, the Sasaki's result is incredible since 
he did not include the exchange diagrams.) Two recent groups indicate 
no evidence for the low-lying pentaquark narrow resonance, 27 ' 36 and one 
recent study suggests a negative-parity pentaquark state in more highly- 
excited region around 1.8GeV. 37 In this section, we perform the accurate 
mass measurement of the 5Q system in anisotropic lattice QCD, and apply 
a new technique to distinguish a compact resonance and a scattering state. 
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2.1. Strategy for High Precession Measurements in Lattice 
QCD 

As a difficulty on the lattice study of the + , several NK scattering states 
appear apart from the compact 5Q resonance state. In this paper, we use 
the term of the + only for the compact 5Q resonance to distinguish it from 
the NK scattering state. In order to examine whether the low-lying 5Q 
state appears as a compact resonance + , we perform the accurate lattice 
QCD calculations with adopting the following three advanced methods. 27 



2.1.1. Usage of Anisotropic Lattice QCD 

We use the anisotropic lattice, where the temporal lattice spacing at is much 
finer than the spatial one a s as shown in Fig. 2. In lattice QCD, hadron 
masses are calculated from the asymptotic temporal behavior of the hadron 
correlators. On the anisotropic lattice, we can get the detailed information 
on the temporal behavior of the 5Q correlator, and hence we can perform 
accurate mass measurements for the low-lying 5Q system. 




Figure 2. Schematic figures of the isotropic lattice (left) and the anisotropic lattice 
(right). On the anisotropic lattice, the temporal lattice spacing at is taken to be smaller 
than the spatial one a s . 



2.1.2. Usage of the Non-N K-type Interpolating Field Operator 

We use a non-N K-type interpolating field to extract the O + (1540) state. 
This choice of the interpolating field would be important and effective. In 
Ref.36, the authors used the NK-type interpolating field and only obtained 
the NK scattering state instead of the compact 5Q state. However, their 
null result may be merely due to a small amount of the compact 5Q com- 
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ponent in the NK -type interpolating field, because their calculation suffers 
from a large contamination of NK scattering states. 
We adopt the non-iVif-type interpolating field, 16 

O a =e abc e ade e b f g (u^C-i 5 d e )(u T } Cd g ){Csl) a , (1) 

for the 5Q state with spin J = 1/2 and isospin 1 = 0. Here a denotes the 
Dirac index, and roman indices a-g are color indices. C = 7472 denotes the 
charge conjugation matrix. Note that the non-iVif-type operator in Eq.(l) 
cannot be decomposed into N and K in the nonrelativistic limit and its 
coupling to the NK state is rather weak. Hence, the 5Q resonance state 
+ can be singled out as much as possible in the present calculation, and 
the results are less biased by the contamination from NK scattering states. 

2.1.3. Application of the Hybrid Boundary Condition Method 

To distinguish compact resonances from scattering states, we propose a 
new method with the "hybrid boundary condition" (HBC) instead of the 
ordinary periodic boundary condition. 27 In the HBC, we impose the anti- 
periodic boundary condition for quarks (u, d) and the periodic boundary 
condition for antiquarks (s), as shown in Table 1. By applying the HBC on 
a finite-volume lattice, the NK threshold is raised up, while the mass of a 
compact 5Q resonance + is almost unchanged. Therefore, we can distin- 
guish a compact 5Q state 0+ from an NK scattering state by comparing 
between the HBC and the standard periodic boundary condition. 

Table 1. The hybrid boundary condition (HBC) to distin- 
guish a compact 5Q resonance ©+ and an NK scattering 
state for the uudds system. The standard boundary condi- 



tion (BC) is alsc 


shown for com] 


>arison. 






u quark 


d quark 


s quark 


HBC 


anti-periodic 


anti-periodic 


periodic 


standard BC 


periodic 


periodic 


periodic 



In lattice QCD with the finite spatial volume L 3 , the spatial momenta 
are quantized as pi = In^ /L (m € Z) under the periodic boundary con- 
dition and pi = (2ni + l)n/L under the anti-periodic boundary condition. 
In the periodic boundary condition, N and K can have zero momenta 
|Pmin| = in the s-wave NK scattering state. The HBC imposes the 
anti-periodic boundary condition for u and d quarks and periodic bound- 
ary condition for s quark, while the periodic boundary condition is usually 
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employed for all u, d, s quarks. In the HBC, the net boundary condi- 
tions of both N(uud, udd) and K(us, ds) are anti-periodic. Then, under 
the HBC, N and K have minimum momenta pi — \/3tt/L in a finite box 
with L 3 , and the threshold of the s-wave NK scattering state is raised up 
as \J m 2 N + p^ in + + Pmm- in contrast to N and K, the compact 5Q 

resonance + (uudds) contains even number of u and d quarks, and hence 
its mass does not shift in the HBC. (See Table 2.) 

Table 2. The net boundary condition for 0+ (uudds), N (uud or 
udd) and K (ds or us) in the hybrid boundary condition (HBC) 
and in the standard boundary condition (BC). 

0+ (uudds) N(uud or udd) K(ds or us) 
HBC periodic anti-periodic anti-periodic 

standard BC periodic periodic periodic 



2.2. Lattice QCD Setup for the Pentaquark Mass 

To generate gluon configurations, we use the standard plaquette action on 
the anisotropic lattice as 27 

^ = |^ E ReTr{l-P ij ( S )} + A 7G ^ ReT r{l-P i4 ( S )},(2) 

s,i<j<3 s,i<3 

with (3 = 2N c /g 2 , the plaquette P a u{s) and the bare anisotropy jq. 

For the quark part, we adopt the O(o)-improved Wilson (clover) fermion 
action on the anisotropic lattice as 

S F = Y,$(x)K(x,vMv), (3) 

with the quark kernel K(x, y) as 
K(x,y) = S x . y - K t {(1 - 74 ) U 4 (x) S x+i y + (1 + 74 ) U\{x - 4) 5 x _ ly ) 
-n s J2{(r- 7i) Ui(x) 6 x+ly + (r + 7i ) u}{x - i) 6 x _ iy } 

i 

— n s ce E ViiGaSx.y — r k s cb E a ijGijS x . y , (4) 

i i<j 

where k s and Kt denote the spatial and temporal hopping parameters, re- 
spectively, denotes the field strength, which is defined through the 
standard clover-leaf-type construction. The Wilson parameter r and the 
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clover coefficients, ce and cb, are fixed by the tadpole-improved tree-level 
values as r = l/£, ce = l/(u s u^) and Cb = where u s and Ut de- 

note the mean-field values of the spatial and the temporal link- variables, 
respectively. 

For the lattice QCD simulation, we use (3 = 5.75 and 12 3 x 96 with the 
renormalized anisotropy a s /a t — 4, which corresponds to 7q = 3.2552. 
In this calculation, the lattice spacing is found to be a s ~ 0.18fm ~ 
(1.1 GeV)- 1 and a t ~ 0.045fm ~ ^GeV)" 1 . We adopt four values of the 
hopping parameter as k = 0.1210 (0.0010) 0.1240 for u and d quarks, and 
use K s -quark = 0.1240 for the s quark. We calculate typical hadron masses 
at each n as shown in Table 3, and find K p h ys . — 0.1261 corresponding to 
the physical situation of m x ~ 0.14 GeV. 



Table 3. Masses of n, p, K and TV for each hopping parameter 
K in the physical unit of GeV. K p h ys . — 0.1261 corresponds to the 
physical situation of ~ 0.14 GeV. 



K 


0.1210 


0.1220 


0.1230 


0.1240 


^phys. 


m w 
m p 
rn K 
m N 


1.005(2) 
1.240(3) 
0.845(2) 
1.878(5) 


0.898(2) 
1.161(3) 
0.785(2) 
1.744(5) 


0.784(2) 
1.085(4) 
0.723(2) 
1.604(5) 


0.656(3) 
1.011(5) 
0.656(3) 
1.460(6) 


0.140 
0.850(7) 
0.530(4) 
1.173(9) 



2.3. Lattice QCD Results for the 0+ (1540) 

Now, using anisotropic lattice QCD, we perform the accurate mass measure- 
ment of the low-lying 5Q states with J = 1/2 and I = in both positive- 
and negative-parity channels from the correlator of the non-iVi-L-type 5Q 
operator with parity projection. 27 

In Fig. 3, we show the lattice QCD results for the masses of lowest 
positive- and negative-parity 5Q states under the standard periodic bound- 
ary condition. After the chiral extrapolation, the lowest positive-parity 5Q 
state is found to be rather heavy as msq^ = \ ) ~ 2.24 GeV, which 
seems to be too heavy to be identified as the + (154O). 

On the other hand, we get a lower mass for the negative-parity 5Q state 
as m 5 Q(J 7r = | ) ~ 1.75 GeV after the chiral extrapolation. This value 
TO5Q ~ 1.75 GeV seems to be closer to the experimental result of m e + — 
1.54 GeV. At this stage, however, this lowest negative-parity 5Q state may 
be merely an NK scattering state, although the non-Aif-type 5Q operator 
used in this calculation includes only a small amount of the NK component. 
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Figure 3. The lowest mass m^Q of the positive- and negative-parity 5Q states plotted 
against m^. The open symbols denote the direct lattice QCD data for positive-parity 
(triangles) and negative-parity (circles). The solid symbols denote the results of the 
chiral extrapolation. The dotted lines indicate the NK thresholds for p-wavc (upper) 
and s-wave (lower) cases. 

To clarify whether the observed low-lying 5Q state is a compact 5Q 
resonance + or an NK scattering state, we apply the new method with 
the hybrid-boundary condition (HBC), and compare the lattice results with 
the HBC and those with the standard periodic boundary condition (BC). 
Recall that, in the HBC, the NK threshold is largely raised up, while the 
mass of the compact 5Q resonance (uudds) is to be almost unchanged, as 
was mentioned in Sect. 2.1.3. 

In Fig. 4, we show the mass of the lowest-lying negative-parity 5Q state 
in lattice QCD with the standard periodic BC and the HBC at each k. 
The symbols denote the lattice QCD results for the 5Q state and the lines 
denote the NK threshold at each k. The left and right figures show the 
results with the standard periodic BC and the hybrid boundary condition 
(HBC), respectively. Note that the NK threshold is estimated to be raised 
up about 200-250 MeV in the HBC. 

As a lattice QCD result, the mass of the 5Q state is largely raised in 
the HBC in accordance with the NK threshold, which indicates that the 
lowest 5Q state observed on the lattice is merely an s-wave NK scattering 
state. In other words, if there exists a compact 5Q resonance + below 
1.75GeV, it should be observed in this lattice calculation with the non-NK- 
type operator, and its mass should be almost unchanged also in the HBC. 
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Figure 4. Comparison between the standard periodic boundary condition (Standard 
BC) and the Hybrid Boundary Condition (HBC) for the lowest mass of the negative- 
parity 5Q system. At each k, the lattice QCD result (the solid symbol) is raised up 
in accordance with the corresponding NK threshold (the solid line). This behavior 
indicates that the low-lying negative-parity 5Q state observed in lattice QCD is an NK 
scattering state rather than a compact 5Q resonance ©+. 



However, there is no such a 5Q state observed in the lattice calculation, 
which means absence of the compact 5Q resonance 6 + below 1.75GeV. 

To conclude, the present lattice QCD calculation at the quenched level 
seems to indicate absence of the low-lying compact 5Q resonance 8 + with 
J = 1/2 and 1 = near 1.54GeV. 27 



2.4. Discussion on Null Result of 0(1540) in Lattice QCD 

In this subsection, we consider the physical consequence of the present 
null result on the low-lying 5Q resonance + in lattice QCD. One possible 
answer is absence of the pentaquark resonance O + (1540), as is indicated 
by recent some high-energy experiments. However, if the 9 + (1540) really 
exists, we deduce the following possibilities. 



2.4.1. Quenching Effects 

First, the present lattice simulation has been done at the quenched level, 
where dynamical quark effects are suppressed. This quenching effect is 
not clear and then it may cause the 5Q resonance G + to be heavier as 
an unknown effect. Actually, recent one lattice study indicates a negative- 
parity 5Q resonance 9 + in a higher-excited region of about 1.8GeV. 37 
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2.4.2. Nontrivial Quantum Numbers of the 0(1540) 

Second, we have considered only the 5Q state with spin J = 1/2 and 
isoscalar 1 = 0. However, the O + (1540) may have other quantum 
numbers, 38,39,40 e.g., spin J = 3/2, isovector I = 1 or isotensor 1 = 2. 
Considering such a possibility, our group started to investigate the 5Q sys- 
tem with higher spin J = 3/2 in lattice QCD as a next step. 

2.4.3. More Complicated Structure of the 6+(1540) 

Third, we have used a localized 5Q interpolating field in this lattice QCD 
calculation. However, the actual + (154O) may have more complicated 
structure beyond the localized interpolating field. Such a possibility has 
been pointed out in the theoretical side. For instance, Karliner and 
Lipkin 41 proposed the diquark-triquark (qq-qqq) picture for the 6, and 
Bicudo et al. 42 pointed out the possibility of the heptaquark picture, where 
the + is described as a bound state of it, K and N. If the 9 + (1540) has 
such a complicated structure, we have to use the corresponding nonlocal 
interpolating field to get its proper information. 

2.5. Necessity of the Wave Function of Pentaquarks 

So far, we have performed the direct mass measurement of 5Q states in 
lattice QCD, where the path integral over arbitrary states is numerically 
calculated on a supercomputer. In the path-integral formalism, however, 
it is rather difficult to extract the state information, such as the wave- 
function of the 5Q state, and therefore only limited simple information can 
be obtained in the direct lattice-QCD calculation. 

Actually, to distinguish the compact 5Q resonance + from NK scat- 
tering states was very difficult in lattice QCD, and hence we had to develop 
a new method with the hybrid boundary condition (HBC). In this respect, 
if the wave function is obtained, one easily finds out whether it is a compact 
resonance state + or not. 

Indeed, to get the wave function is very important to clarify the further 
various properties of the 5Q state such as the underlying structure and 
the decay width, which cannot be obtained practically only with the direct 
lattice-QCD calculation. 

Then, apart from the direct lattice-QCD calculation, we have to seek the 
way to obtain the proper wave function of the 5Q state. To do so, we need a 
proper Hamiltonian for the 5Q system based on QCD. One possible way in 
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this direction is to construct the quark model from QCD, as was mentioned 
in Section 1. In the next section, we study the inter-quark interaction in 
multi-quark systems directly from QCD, and aim to construct the QCD- 
based quark-model Hamiltonian. 

3. Lattice QCD Study for the Inter-Quark Interaction in 
Multi-Quark Systems 

In this section, we perform the first study of the inter-quark interaction in 
multi-quark systems in lattice QCD, 22,23,24 ' 25,26 and seek for the QCD- 
based quark-model Hamiltonian to describe multi-quark hadrons. The 
quark-model Hamiltonian consists of the kinetic term and the potential 
term, which is not known form QCD in multi-quark systems. 

As for the potential at short distances, the perturbative one-gluon- 
exchange (OGE) potential would be appropriate, due to the asymptotic 
nature of QCD. For the long-range part, however, there appears the con- 
finement potential as a typical non-perturbative property of QCD, and its 
form is highly nontrivial in the multi-quark system. 

In fact, to clarify the confinement force in multi-quark systems is one 
of the essential points for the construction of the QCD-based quark-model 
Hamiltonian. Then, in this paper, we investigate the multi-quark potential 
in lattice QCD, with paying attention to the confinement force in multi- 
quark hadrons. 

3.1. The Three-Quark Potential in Lattice QCD 

So far, only for the simplest case of static QQ systems, detailed lattice 
QCD studies have been done, and the QQ potential Vqq is known to be 
well described by the Coulomb plus linear potential as 17 > 18 > 19 > 20 

^QqM— ^+<W + C QQ (5) 

with r being the inter-quark distance. 

To begin with, we study three-quark (3Q) systems in lattice QCD to 
understand the structure of baryons at the quark-gluon level. Similar to 
the derivation of the QQ potential from the Wilson loop, we calculate the 
3Q potential V^q from the 3Q Wilson loop Wzq in SU(3) lattice QCD with 
(/3=5.7, 12 3 x 24), 08=5.8, 16 3 x 32), 08=6.0, 16 3 x 32) and (/? = 6.2, 24 4 ) 
at the quenched level. For more than 300 different patterns of spatially- 
fixed 3Q systems, we perform accurate and detailed calculations for the 3Q 
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potential, 19 ' 20 ' 21 ' 23 ' 25 and find that the lattice QCD data of the 3Q potential 
Vjjq are well described by the Coulomb plus Y-type linear potential, i.e., 
Y-Ansatz, 

V 3Q = -A 3Q £ r— — I + KqL 3 ^ + C 3Q: (6) 

, ■ 1 i 1 j\ 

within 1%-level deviation. 19 ' 20 ' 21 ' 23 ' 25 Here, L^ n is the minimal total length 
of the color flux tube, which is Y-shaped for the 3Q system. 

To demonstrate the validity of the Y-Ansatz, we show in Fig. 5 the lattice 
QCD data of the 3Q confinement potential V 3 C Q nf , i.e., 3Q potential sub- 
tracted by the Coulomb part, plotted against Y-shaped flux-tube length 
-^min- F° r each p, clear linear correspondence is found between 3Q confine- 
ment potential V 3 Q nf and L^ n , which indicates the Y-Ansatz for the 3Q 
potential. 23 ' 25 
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Figure 5. The lattice QCD result for the 3Q confinement potential Vjjq , i.e., the 3Q 
potential subtracted by its Coulomb part, plotted against Y-shaped flux-tube length 
^mhi a * 0=5-8, 6.0 and 6.2 in the lattice unit. The clear linear correspondence between 
3Q confinement potential V r 3 C Q nf and indicates the Y-Ansatz for the 3Q potential. 



Here, we consider the physical meaning of the Y-Ansatz. Apart from an 
irrelevant constant, the Y-Ansatz, Eq.(6), consists of the Coulomb term and 
the Y-type linear potential, which play the dominant role at short and long 
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distances, respectively. The Coulomb term would originate from the one- 
gluon-exchange (OGE) process. In fact, at short distances, perturbative 
QCD is applicable, and therefore the inter-quark potential is expressed as 
the sum of the two-body one-gluon-exchange (OGE) Coulomb potential. 
The appearance of the Y-type linear potential supports the flux-tube 



picture 



43,44.45.46 



at long distances, where there appears the color flux tube 



linking quarks inside hadrons with its length minimized. In particular, the 
confinement force in baryons corresponds to the Y-shaped flux tube, which 
implies existence of the three-body interaction in baryons. 

In usual many-body systems, the main interaction is described by a two- 
body interaction and the three-body interaction is a higher-order contribu- 
tion. In contrast, as is clarified by our lattice-QCD study, the quark con- 
finement force in baryons is a genuinely three-body interaction, 19,20 which 
is one of significant features of QCD. In fact, the appearance of the Y- 
type junction and the three-body confinement force reflect the SU(3) group 
structure in QCD, e.g., the number of color, N c — 3, and is peculiar to 
QCD. 19,20 In this sense, the study of the 3Q system is very important to 
get a deeper insight of the QCD physics. 




Figure 6. Lattice QCD results for Y-shaped flux-tube formation in the 3Q system: a 
bird's-eye view and a contour map of the action density of the spatially-fixed 3Q system 
in maximally-Abelian projected QCD. The distance between the junction and each quark 
is about 0.5 fm. 



Recently, a clear Y-type flux-tube formation is actually observed for 
spatially-fixed 3Q systems in lattice QCD as shown in Fig. 6. 23 ' 25, 47 Thus, 
together with recent several other analytical and numerical studies, 48 ' 49,50 
the Y-Ansatz for the static 3Q potential is almost settled. This result 
indicates the color-flux-tube picture for baryons. 
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3.2. The OGE Coulomb plus Multi-Y Ansatz 

Now, we proceed to multi-quark systems. We first consider the theoretical 
form of the multi-quark potential, since we will have to analyze the lattice 
QCD data by comparing them with some theoretical Ansatz. 

By generalizing the lattice QCD result of the Y-Ansatz for the three- 
quark potential, we propose the one-gluon-exchange (OGE) Coulomb plus 
multi-Y Ansatz, 22 > 23 > 24 > 25 < 26 

2 x a T a 

i<j 11 -? 1 

for the potential form of the multi-quark system. Here, the confinement 
potential is proportional to the minimal total length i m ; n of the color flux 
tube linking the quarks, which is multi-Y shaped in most cases. 

In the following, we study the inter-quark interaction in multi-quark 
systems in lattice QCD, and compare the lattice QCD data with the theo- 
retical form in Eq.(7). Note here that the lattice QCD data are meaningful 
as primary data on the multi-quark system directly based on QCD, and do 
not depend on any theoretical Ansatz. 



3.3. Formalism of the Multi-Quark Wilson Loop 

Next, we formulate the multi-quark Wilson loop to obtain the multi-quark 
potential in lattice QCD. 22 ' 23 ' 24 ' 25 ' 26 

Similar to the derivation of the QQ potential from the Wilson loop, the 
static multi-quark potential can be derived from the corresponding multi- 
quark Wilson loop. We construct the tetraquark Wilson loop W 4 q and 
the pentaquark Wilson loop W5Q in a gauge invariant manner as shown in 
Figs. 7(a) and (b), respectively. 

The tetraquark Wilson loop W 4 q and the pentaquark Wilson loop W5Q 
are defined by 

W 4Q = -tr(Mi^i 2 M 2 £i2), 

W 5Q = l-e abc e a ' b ' c ' M aa ' (R 3 R 12 R 4 ) bb ' (L 3 L 12 L 4 ) CC ' , (8) 
where M, Mj, Lj and Rj (i=l,2, j=l,2,3,4) are given by 

M, Mi, Rj,Lj = Pcxp{ig I dx^A^x)} e SU(3) C . (9) 
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Figure 7. (a) The tetraquark Wilson loop W4Q. (b) The pentaquark Wilson loop W5Q. 
The contours M, Mj, Rj, Lj(i = 1,2, j = 3,4) are line-like and Rj,Lj(j = 1,2) are 
staple-like. The multi-quark Wilson loop physically means that a gauge-invariant multi- 
quark state is generated at t = and annihilated at t = T with quarks being spatially 
fixed in R 3 for < t < T. 



Here, R12 and ii 2 are denned by 



j~ta a — 1 ^abc^a'b'c' r?bb' tjcc' 



a'a — ^abc^a'b' c' r bb' r cc' 



La a — ^ 
12 = ~ e 



j t c 



(10) 



The multi-quark Wilson loop physically means that a gauge- invariant multi- 
quark state is generated at t — and annihilated at t = T with quarks being 
spatially fixed in R 3 for < t < T. 

The multi-quark potential is obtained from the vacuum expectation 
value of the multi-quark Wilson loop as 



V 4Q = - lim ^In<W4Q), V 5Q = - km ^ln(W 5Q ). (11) 

1 — *oo 1 1 — >oo 1 



3.4. Lattice QCD Setup for the Multi-quark Potential 

Here, we briefly summarize the lattice QCD setup in this calculation. For 
the study of the multi-quark potential, the SU(3) lattice QCD simulation 
is done with the standard plaquette action at (3 = 6.0 on the 16 3 x 32 
lattice at the quenched level. (The calculation for large-size multi-quark 
configurations are performed by identifying 16 2 x 32 as the spatial size.) 

In this calculation, the lattice spacing a is estimated as a ~ 0.104fm, 
which leads to the string tension CTqq = (427MeV) 2 in the QQ potential. 23 
We use 150 gauge configurations for the 5Q potential simulation and 300 
gauge configurations for the 4Q potential simulation. The smearing method 
is used for the enhancement of the ground-state component. We here adopt 
a = 2.3 and the iteration number iV smr = 40, which lead to a large enhance- 
ment of the ground-state component. 22 > 23 > 24 > 25 > 26 
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3.5. Lattice QCD Result of the Pentaquark Potential 

We perform the first study of the pentaquark potential V5Q in lattice QCD 
for 56 different patterns of QQ-Q-QQ type pentaquark configurations, as 
shown in Fig. 8. As the conclusion, the lattice QCD data of V5Q are found 
to be well described by the OGE Coulomb plus multi-Y Ansatz, i.e., the 
sum of the OGE Coulomb term and the multi-Y-type linear term based on 
the flux-tube picture. 22,23,24,25 



3 3 




Figure 8. A QQ-Q-QQ type pentaquark configuration. In the 5Q system, (Qi, Q2) and 
(Q31Q4) form 3 representation of SU(3) color, respectively. The lattice QCD results 
indicate the multi- Y-shaped flux-tube formation in the QQ-Q-QQ system. 



We show in Fig. 9 the lattice QCD results of the 5Q potential V5Q for 
symmetric planar 5Q configurations as shown in Fig. 8, where each 5Q sys- 
tem is labeled by d = Q1Q2/2 and h = Q1Q3. 




h Lmin 



Figure 9. Lattice QCD results of the pentaquark potential Vsq for symmetric planar 
5Q configurations in the lattice unit: (a) V 5 q v.s. (d, h) and (b) V5Q v.s. L^ n . The 
symbols denote the lattice QCD data, and the curves the theoretical form of the OGE 
plus multi-Y Ansatz. 



In Fig. 9, we add the theoretical curves of the OGE Coulomb plus multi-Y 
Ansatz, where the coefficients (A$q, 05Q) are set to be (A3Q, 0-3 q) obtained 
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from the 3Q potential. 20 (Note that there is no adjustable parameter in the 
theoretical Ansatz apart from an irrelevant constant.) In Fig. 9, one finds a 
good agreement between the lattice QCD data of V5Q and the theoretical 
curves of the OGE Coulomb plus multi-Y Ansatz. 

In this way, the pentaquark potential V5Q is found to be well described 
by the OGE Coulomb plus multi-Y Ansatz as 22,23,24,25 

^5Q = - A 5Q {(— + —) + - — + — + — + — 

t-12 r 34 2 n 5 r 25 r 35 r 45 

+ \(— + — + — + — )} + tr BQ L^ n + C BQ , (12) 
4 n 3 n 4 r 23 r 24 

where is the distance between Q, and Qj. Here, L^ n is the minimal 
total length of the flux tube, which is multi- Y-shaped as shown in Fig. 8. 
This lattice result supports the flux-tube picture for the 5Q system. 

3.6. Tetraquark Potential and Flip-Flop in Lattice QCD 

We study the tetraquark potential V4Q in lattice QCD for about 200 differ- 
ent patterns of QQ-QQ configurations, and find the following results. 23 > 24 > 26 

(1) When QQ and QQ are well separated, the 4Q potential V4Q is well 
described by the OGE Coulomb plus multi-Y Ansatz, which indi- 
cates the multi- Y-shaped flux-tube formation as shown in Fig. 10(a). 

(2) When the nearest quark and antiquark pair is spatially close, the 
4Q potential V4Q is well described by the sum of two QQ potentials, 
which indicates a "two-meson" state as shown in Fig. 10(b). 

Qi Q 3 
• o 



Q 4 



Figure 10. (a) A connected tetraquark (QQ-QQ) configuration and (b) A "two-meson" 
configuration. The lattice QCD results indicate the multi- Y-shaped flux-tube formation 
for the connected 4Q system. 

We show in Fig. 11 the lattice QCD results of the 4Q potential V4Q for 
symmetric planar 4Q configurations as shown in Fig. 10, where each 4Q 
system is labeled by d = Q1Q2/2 and h = Q1Q3. 
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Figure 11. Lattice QCD results of the tetraquark potential V4Q for symmetric planar 
4Q configurations in the lattice unit. The symbols denote the lattice QCD data. The 
solid curve denotes the OGE plus multi-Y Ansatz, and the dotted-dashed curve the 
two-meson Ansatz. 

For large value of h compared with d, the lattice data seem to coincide 
with the solid curve of the OGE Coulomb plus multi-Y Ansatz, 



(13) 
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^4Q — -AlQ 

+ a iQ L tSn + C4Q, 

where L^ n is the minimal total length of the flux tube, which is multi-Y- 
shaped as shown in Fig. 10(a). Here, the coefficients (A4q,04q) are set to 
be (A 3 q,(73q) obtained from the 3Q potential. 20 

For small h, the lattice data tend to agree with the dotted-dashed curve 
of the "two- meson" Ansatz, where the 4Q potential is described by the sum 
of two QQ potentials as VQq(ri 3 ) + V Q Q(r 2 4) = 2Vqq(/j). 

Thus, the tetraquark potential V4Q is found to take the smaller energy 
of the connected 4Q state or the two-meson state. In other words, we 
observe a clear lattice QCD evidence of the "flip-flop", i.e., the flux-tube 
recombination between the connected 4Q state and the two-meson state. 
This lattice result also supports the flux-tube picture for the 4Q system. 



3.7. Proper Quark-Model Hamiltonian for Multi-Quarks 

From a series of our lattice QCD studies 19 ' 20 ' 21 ' 22 ' 23 ' 24 ' 25,26 on the inter- 
quark potentials, the inter-quark potential is clarified to consist of the one- 
gluon-exchange (OGE) Coulomb part and the flux-tube-type linear confine- 
ment part in QQ-mesons, 3Q-baryons and multi-quark (4Q, 5Q) hadrons. 

Furthermore, from the comparison among the QQ, 3Q, 4Q and 5Q 
potentials in lattice QCD, we find the universality of the string tension a, 

a QQ — a 3Q — a 4:Q — 05Q, (14) 
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and the OGE result of the Coulomb coefficient A as 

2^QQ - ^3Q - Mq ^ ^5Q 
in Eqs.(5), (6), (12) and (13). 



(15) 



Here, the OGE Coulomb term is considered to originate from the OGE 
process, which plays the dominant role at short distances, where pertur- 
bative QCD is applicable. The flux-tube-type linear confinement would 
be physically interpreted by the flux-tube picture, where quarks and anti- 
quarks are linked by the one-dimensional squeezed color-electric flux tube 
with the string tension a. 

To conclude, the inter-quark interaction would be generally described 
by the sum of the short-distance two-body OGE part and the long-distance 
flux-tube-type linear confinement part with the universal string tension 
cr ~ 0.89 GeV/fm. 

Thus, based on the lattice QCD results, we propose the proper quark- 
model Hamiltonian H for multi-quark hadrons as 



where L m ; n is the minimal total length of the flux tube linking quarks. Vq GE 
denotes the OGE potential between ith and jth quarks, which becomes 
the OGE Coulomb potential in Eq.(7) in the static case. M, denotes the 
constituent quark mass. The semi-relativistic treatment would be necessary 
for light quark systems. 

It is desired to investigate various properties of multi-quark hadrons 
with this QCD-based quark model Hamiltonian H. 

4. Summary and Concluding Remarks 

We have studied tetraquark and pentaquark systems in lattice QCD, moti- 
vated by the recent experimental discoveries of multi-quark candidates such 
as the 6+(1540). 

First, we have performed accurate mass calculations of low-lying 5Q 
states with J = 1/2 and I = in both positive- and negative-parity chan- 
nels in anisotropic lattice QCD. We have found that the lowest positive- 
parity 5Q state has a large mass of about 2.24GeV after the chiral extrap- 
olation. To single out the compact 5Q state from NK scattering states, we 
have developed a new method with the hybrid-boundary condition (HBC), 
and have found no evidence of the compact 5Q state below 1.75GeV in the 
negative-parity channel. 




(16) 
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Second, we have performed the first study of the multi-quark potential 
in lattice QCD to clarify the inter-quark interaction in multi-quark sys- 
tems. We have found that the 5Q potential V$q for the QQ-Q-QQ system 
is well described by the "OGE Coulomb plus multi-Y Ansatz" : the sum of 
the one-gluon-exchange (OGE) Coulomb term and the multi-Y-type linear 
term based on the flux-tube picture. The 4Q potential V4Q for the QQ- 
QQ system is also described by the OGE Coulomb plus multi-Y Ansatz, 
when QQ and QQ are well separated. On the other hand, the 4Q system is 
described as a "two-meson" state with disconnected flux tubes, when the 
nearest quark and antiquark pair is spatially close. We have observed a 
lattice-QCD evidence for the "flip-flop", i.e., the flux-tube recombination 
between the connected 4Q state and the "two-meson" state. On the confine- 
ment mechanism, we have clarified the flux-tube-type linear confinement in 
multi-quark hadrons. Finally, we have proposed the proper quark-model 
Hamiltonian based on the lattice QCD results. 

Acknowledgements: F.O. thanks the organizers of QNP05 for their 
warm hospitality in Korea. The lattice QCD Monte Carlo calculations have 
been performed on NEC-SX5 at Osaka University and SR8000 at KEK. 
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